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§ 1 Introduction 

Lieb[I]l3] and Sutherland [3] solved the six- vertex model which is a lattice 
model now known to be equivalent to the XXZ quantum spin chain. The 
six-vertex model satisfies the Yang-Baxter equation which plays fundamental 
role in integrable lattice models. [HI E] 

Halsey et al. [7] analyzed the spectrum of singularities lying upon possibly 
fractal sets, for example those found in dynamical systems theory. This 
was generally formulated by Edgar and Mauldin[8j as a theory of new sub- 
structure lying in fractal sets. This sub-structure is called " multifractal" . 

In this paper, the fractal structure lying in the six-vertex model is for- 
mulated as a fractal generated by an iterated function system (IFS). We 
find that the notions which correspond to the transfer matrix theory and 
an equivalence relation of boundary conditions in lattice systems have also 
been introduced in the area of fractal geometry. The functional relation 
for the fractal dimension of the fractal set corresponding to the six-vertex 
model shows structural change in the thermodynamic limit, depending on an 
anisotropy parameter. 

This correspondence can be formulated for systems suitable to the transfer 
matrix treatment. This gives us a possibility to find relations between theo- 
ries of fractal sets and those of lattice models. In particular, models which 
satisfy the Yang-Baxter relation provide fractal sets that can be mapped to 
integrable systems. This is especially interesting when we consider the fact 
that there exist infinite hierarchy of solvable models [H [10] though we do 
not have sufficiently many examples of fractal sets where we can obtain the 
fractal dimensions exactly. The relation is also interesting from the point of 
view that the six- vertex model is governed by the quantum group. [B] The 
correspondence itself is also valid if one considers nonintegrable cases. 



This note is arranged as follows. |§ 2 [ is a short review of the six- vertex 
model and the equivalence relation that classifies the free energy. § 3 is 
also a short review of the IFS. Informed readers in respective fields may skip 
these sections. The correspondence between the transfer matrix method and 
IFS, the ra-equivalence and the transitivity condition of IFS fractals, and the 
relation between the fractal dimension and the free energy of the six-vertex 



model is explained in § 4 
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2 The six- vertex model and the n-equivalence 
relation 



Let us introduce the six- vertex model. We consider a rectangular lattice with 
h rows and w columns, and assign an arrow to each bond with the rule that 
two of four arrows point in and other two point out at each site (Fig. 1(a)). 
Then six types of local arrow arrangements are possible. The site together 
with the four arrows around it is called the vertex. Each vertex is assumed 
to have finite and non- negative energies ei, €2 or €3. The vertex energies are 
assumed to be invariant under the reversal of all arrows. 

When a line is drawn on each bond that points down or points left, 
a one-to-one correspondence between the configurations of arrows and the 
configurations of lines on the lattice can be found (Fig. 1(a) and (b)). These 
lines do not intersect each other. Each line begins from one bond on the 
boundary, and continues until it reaches another bond on the boundary. 
Thus the number of lines on the lattice is determined by the number of lines 
on the boundary. 

The free energy of the six-vertex model on the rectangle was first ob- 
tained with the cyclic boundary conditions in two directions. [UH] The trans- 
fer matrix V of the six-vertex model and the Hamiltonian Tixxz of the XXZ 
quantum spin chain commute, and hence they share the same eigenvectors. 
There also exists a simple relation between the eigenvalues of V and those 
of Ti-xxz (see for example [6]). 

Next let us introduce an equivalence relation of boundary conditions 
which is called the n-equivalence. |Tl] Let us consider models in which each 
of the local variables takes one of a finite number of discrete states. We con- 
sider the lattice where the number of boundary sites A^' is of a lower order 
than the total number of sites N in the thermodynamic limit: A^' = o{N) 
{N 00). 

Let us consider a site on the lattice and the number of steps (the number 
of bonds) n' that is necessary to reach one of the sites on the boundary. 
There exists the minimum of n' for each site. Then let us consider the sites 
where the minimum of n' is equal to n. We call them the ra-boundary sites. 
Let us consider the set of bonds between {n — 1)- and ra-boundary sites, and 
call them the ra-boundary bonds. The set of n-boundary sites together with 
the n-boundary bonds is called the ra-boundary, and configurations on the 
n-boundary are called the n-boundary configurations. 
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Let us introduce a set of n-boundary configurations {Fj} which is the set 
of all the allowed configurations on the n-boundary under a specific boundary 
condition F on the actual boundary of the lattice. Then the equivalence of 
boundary conditions is defined as follows: boundary conditions F and F' are 
n-equivalent if {F,} = {F^} as a set of ra-boundary conditions. 

In some models such as the six-vertex model, the free energy depend 
on the boundary condition still in the thermodynamic limit. However it is 
derived p!T| [12j that the free energies with boundary conditions F and F' are 
identical in their thermodynamic limit, if the boundary conditions F and F' 
are n-equivalent with a finite n throughout the limit. It is also true when 
the number n diverges but satisfies n = o{N/N') {N — * oo). 

With the use of this equivalence, it is derived [12] that the free energy of 
the six-vertex model on domain D, with fixed density of lines pi for the hori- 
zontal bonds and fixed density of lines p2 for the vertical bonds on the bound- 
ary, are identical to each other in the thermodynamic limit: / = /(pi,P2)- 

The equivalence can be introduced on a part of the boundary, for example 
on the first row of the rectangle. In this paper, we will introduce this type 
of n-equivalence. In this case the corresponding n-boundary sites are those 
on the (n -|- l)-th row of the lattice. 

The n-equivalence corresponds to the irreducibility of the transfer matrix, 
and corresponds to the regularity of the stochastic matrix. It will be noted 
later that the n-equivalence also corresponds to the transitivity condition in 
fractal geometry. 

§ 3 Iterated Function System (IFS) 

Next let us introduce the iterated function system (IFS), which provides a 
way to construct fractal sets through iterations of contractions defined by a 
set of functions. First we will consider the Cantor set which is one of the 
simplest fractal sets generated by an IFS. Let us consider the interval Jq = 
[0, 1] and introduce two contractions defined by the functions Fi{x) = |a; and 
F2{x) = I + ^x. The functions Fi and F2 generate two similar subsets with 
the contraction ratio r = 1/3: Fi([0, 1]) = [0,1/3] and F2([0,l]) = [2/3,1]. 
Thus one obtains h = i^i(/o) U Faj/o) = [0, 1/3] U [2/3, 1]. Beginning from 
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Ii and again operating Fi and F2, one obtains 

h = Fi(/i)UF2(/i) (1) 
= [0, 1/9] U [2/9, 1/3] U [2/3, 7/9] U [8/9,1]. (2) 

Operating the functions Fi and F2 iteratively, one obtains the subset /„ 

4 = Fi(/„_i)UF2(/„_i). (3) 

Taking the hmit n ^ 00, it is known that there remains a non- vanishing 
subset, which we call the Cantor set. The Cantor set has been generated by 
the set of functions {^1,^2}. This procedure can be generalized to the case 
of the set of finite number of functions {Fi, F2, . . . , Fp}, which we call the 
iterated function system. 

One can also assume that an additional index j is assigned to each in- 
terval. Let us introduce contraction functions Fij, which operates only on 
the intervals of type j and generates the intervals of type i. Not all the 
intervals can be generated from the interval of type j: some of the genera- 
tions j ^ i may be prohibited. This restriction is often displayed through 
the graph (Fig. 2). Each arrow from i to j is usually written as a symbol 
for the function Fij. Two or more arrows from i to j, or from i to i, with 
different contraction ratios may exist. This kind of restricted IFS is called 
the graph-directed IFS. [T3lfT5] 

In the case of the Cantor set, 2 similar small intervals are generated from 
an interval in each iteration with the contraction ratio r = 1/3. In the 
case of the simple equal division of a d-dimensional interval, the number of 
generated small intervals should be (1/r)'^. Then the similarity dimension ds 
is introduced through the relation 

2 = (l/r)'^^ (4) 

and one obtains = log 2/ log 3 = 0.6309 ■ ■ ■ < 1. The similarity dimension 
ds, which is introduced in the case of self-similar sets, is the simplest example 
of fractal dimensions. 

One can introduce another fractal dimension. Let E he a non-empty and 
bounded subset of R''. Let Ns{E) be the smallest number of d-dimensional 
intervals (d-dimensional boxes) of diameter 6, with which one can cover the 
set E. Then the box-counting dimension of E is defined as 

dimBi5 = limi^^^iMf). (5) 
.5^0 -log 5 
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This definition means that the smallest number of boxes to cover E is of 
order {1/6Y^, where = dimB-E, in the limit 6—^0. 

The Hausdorff dimension of E, usually written as = dimn-E, is defined 
by introducing a countable collection of open sets to cover E, instead of the 
set of boxes. Taking the limit 5 — > 0, where 6 is the suprimum of the 
diameters of the open sets, we can introduce the Hausdorff measure of E, 
and the Hausdorff dimension is defined as the dimension where the measure 
jumps from cx) to (see details for example in [IS]). The Hausdorff dimension 
might be the most sophisticated dimension to measure fractal. 

In the case of complicated fractal sets, fractal dimensions often take dif- 
ferent values from each other. However in our case, fractal sets generated by 
IFS, du is equal to d^. In the case of self-similar sets, we have d^ = ds = ds- 

§ 4 Fractal structure of the six- vertex model 

Let us consider the six-vertex model on a rectangle with w columns and h 
rows. When we fix the line configuration on the first row of the rectangle, we 
see that not all of the configurations are possible on the next low, because 
of the six- vertex restriction (Fig. 1(b)). The set of allowed configurations 
are determined by the configuration on the first low. Let us assume that the 
configuration is type j on the first low, and that the configurations ii, ^2, • ■ ■ , 

are allowed on the second low. Then assuming the configuration ii on the 
second low, configurations i'l, i'2, ■ ■ ■ , i'p^ are allowed on the third low. The set 
of allowed configurations {i'i,i2, • • • ^'^pg} is determined by the configuration 
ii. The configurations are generated by operating the row to row transfer 
matrix V. This iteration procedure corresponds to the generation of fractal 
sets by means of the graph-directed IFS. This is the simple summary of the 
procedure to relate lattice models and IFS fractals. 

Now let us strictly define the graph-directed IFS and a possible set of 
fractals corresponding to the six-vertex model. 

Let us consider a finite number of "dots" labeled by index j. Let us 
introduce a set of directed edges, where each edge e-^^ starts a dot i and ends 
at a dot j. A pair of dots i and j may be joined by several edges distinguished 
by the index k. Edges from dot i to i itself may also exist. Let us introduce 
a contraction function E^^^: H'^ — > R*^, corresponding to each e[j\ Let r-^^ 

(k) 

be the contraction ratio of F, •;• , which is the infimum of the number r that 
satisfies \E^^\x) — El^^\y)\ < r\x — y\ for all x,y G R". It is assumed that 
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< ri^^ < 1. 

Let us consider a set of n' directed edges (cjfc^, e^j^j, . . . , ek^,_j), ki ^ i,j 

(/ = 1, ... ,n' — 1), which form a sequential path from i to j. Let £>j be 
the set of such sequential n' edges. We assume the transitivity condition, 
i.e. there is a positive integer n which satisfies that, for all i, j, there exists 

(n - ■) 

an integer riij < n such that * is not empty. The transitivity condition 
means that there exist finite sequential paths in the graph joining every pair 
of dots i and j. 

Then it is known that there exists a unique family of non-empty fractal 
sets {Ej} such that 

j k 

The set of functions {Pj^^''} is called a graph-directed IFS, and the fractal 
sets {Ej} are called a family of graph-directed sets. 

When the right-hand side of ([6]) is disjoint, we say that the set of functions 
{Fj^j^} satisfies the (strong) separation condition. 

In the case of the six- vertex model, each "dot" labeled j corresponds to 
the set of line configurations allowed under the condition that the boundary 
line configuration on the first row is fixed and labeled j. Each directed edge 
corresponds to an allowed generation of a configuration, from a configuration 
j on a row of vertical bonds to a configuration i on the next row. 

The functions {-F//'*} are introduced as follows. Let us introduce an order 
of bonds as shown in Fig. 3 and assign numbers {si}, where si = +1 (or 
si = 0) if the arrow on the Z-th bond points down/left (or up/right). When 
we assume w columns and h rows and the cyclic boundary condition in the 
horizontal direction, the number of bonds is equal to 2wh + w. 

The set of numbers {si, . . . , s^} is determined by the line configuration on 
the first w vertical bonds, which is the boundary condition on the first row. 
Let us introduce the number j = Ym=i which works as an index for the 
configuration on the first row of vertical bonds. The correspondence from 
{si, . . . , s^} to + 2'"^] provides a mapping from an allowed configuration 
to an interval. 

When we consider the next 2w bonds, we will find a row of w vertices. 
Let us introduce contractions of each interval [j, j + 2"^"] as 

[j,j+2-]H^U[x2,a;2 + 2-3-r~,,], (7) 

X2 
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where X2 = j + Yl^=w+i sfi^^ provides an index of the hne configuration on 
the lattice up to the second row of vertical bonds, is an weight which 
satisfies < Vij < 1 and will be fixed later in f|T5|l . and in the union |J^^ the 
index X2 runs over all the allowed line configurations on the added 2w bonds 
with fixed {si, S2, ■ ■ ■ , s^}- Let us introduce a function 

hw 

x = J2si2'^^p{x,h)= sa~^''^^-'^-\ (8) 

1=1 l={h-l)w+l 

which is an index that distinguishes line configuration {s(fc_i)^+i, S(^k-i)w+2^ ■ ■ ■ , s^w}, 
e.g. p{x2, 1) = j. The contraction functions F^p are introduced as 

2-) = (2,X2,2-=^-rif), (9) 

where p{x2, 2) = k and p{x2, 3) = i, and the contraction ratio of F^^ is r^j^ = 
2-2«)~{fc)_ Eg^f.]^ index k indicates an allowed configuration of {s^+i, . . . , S2w} 
and the index i indicates an allowed configuration of {s2w+i, S2w+2, ■ ■ ■ , s^w}- 
Thus a finite set of contraction functions is obtained corresponding to all 
the allowed generations from the configuration j to X2 where p{x2, 3) = i. 
After the contractions by F^^^ we find a set of intervals, where each interval 

[x2, X2 + 2~^^rlj^] corresponds to an allowed line configuration {si, . . . , 53^,} 
on the lattice. 

When the number of rows is generally increased as /i ^— /i + 1, 2w bonds 
will be added. The contraction function works as 

F^P{h,x,,A) = (/^ + l,x,+i,2-2-rt)A), (10) 

where Xh+i = Xh + 2~^'^^-^^'"k + 2-^'^"'^, p{xh, 2h ~ I) = p{xh+i, 2h - I) = j, 
p{xh+i,2h) = k and p{xh+i,2h + 1) = i. This provides contractions of 
intervals 

[xh, Xh + A]^ [xH+u Xh+i + 2-2"'^;^ A]. (11) 

Contraction functions F^^^ corresponding to the generation h ^ h + ho {Hq G 
N) can be introduced by straightforward generalizations. 

Taking the limit h 00, we obtain fractal sets {Ej} generated by the 
graph-directed IFS where the set of contraction functions is given by {F^^^}. 
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As for the fractal dimension of graph-directed sets generated by a graph- 
directed IFS {Pj^Jr}, the following results are generally obtained [121 HH and 
see for example [151 [E] ■ Let us introduce a g-dimensional matrix A^^'' 

{A^'\.' = ' ^ (12) 

7 

where r^J^, is the contraction ratio of F^^J, and a,a' = 1,2, ... ,q. Let Pq(s) be 
the eigenvalues of A^'^\ From the Perron- Frobenius theorem, we know that 
there exists a real and positive eigenvalue Pmax(s) which satisfies |pa(s)| < 
Pmax(s) for all a. Then, 

Theorem 1 Let Ei, E2,. ■ ■ , Eg be a family of graph- directed sets generated 
by a graph- directed IFS {Fj^J,} that satisfies the transitivity and the separation 
conditions. Then, 

1) dimn-Ea = dimB-E„ (a = 1, 2, . . . , g), 

2) there exists an positive number s such that 

dimn^^i = dimH£^2 = ■ ■ ■ = dimn-Eg = s, (13) 

3) the number s is the unique solution of the equation 

Pmax(s) = 1. (14) 
In the case of the six-vertex model, let 

= exp(-A)6S;)), (15) 

(k) 

where /3o is a non- negative constant and e\j is the sum of the energy of 
added w vertices. The contraction ratio is r 5f = (l/2)2«'fj^^\ and thus the 
matrix A^^^ becomes A'^^'^ = (1/2)^'"*1^', where V is the row to row transfer 
matrix with the restricted q basis. Because of the first factor (1/2)^"'*, the 
separation condition is satisfied. 

If the configuration j and i are n-equivalent, the configuration i can be 
generated from j within 2n times operations of the transfer matrix. [12j Thus 
the transitivity condition corresponds to the fact that all the boundary con- 
figurations on the first row are ra-equivalent to each other with some finite 
n: when we consider the set of all the n-equivalent boundary conditions, the 
corresponding graph-directed IFS is transitive. 
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When we introduce all the admissible line configurations with fixed m 
lines on the first row, these configurations form the complete set of n-equivalent 
boundary conditions. The transfer matrix is block diagonalized according to 
m. The fractal dimension s is obtained from the condition (flU), i.e. 



i^f'^'K.AsPo) = 1, (16) 

where Aniax(/?) is the maximum eigenvalue of the block element of the transfer 
matrix, with the temperature /3 = s/3o- Hence we obtain the relation 

s = ^^T^, -PfiP) = - log A^a.(/3), (17) 
2 log 2 w 

where /(/?) is the free energy of the six- vertex model with m lines. We already 
knowpT| [T2] that n-equivalent (i.e. transitive) boundary conditions yield 
the identical /(/?) (i.e. identical fractal dimension) in the thermodynamic 
limit, that is consistent with (fT3|) . Eq. (fT7|) is the functional relation for 
s. The dimension s is not identical with the multifractal dimension. The 
correspondences of the concepts in the theory of lattice models and those in 
the fractal geometry are summerised in Table 1. 

The resulted s depends on the 'anisotropy' parameter A of the six-vertex 
model only through the ratio rjj, and hence the correspondence is valid for 
all A. 

However Amax has different functional form depending on A. The maxi- 
mum eigenvalue Amax of the transfer matrix for the six-vertex model is (see 
for example [20j ) 

w w 

Amax(/3) = n ^(^^•) +^"11 ^(^^•)' (IS) 

j=l j=l 

where L{z) = {ab+ {c'^ -b^)z)/{a'^ - abz), M{z) = {a'^ - - abz)/{ab-b'^z), 
a = e~^^^, b = e~^'^'^ and c = e~^^^. The complex numbers Zj satisfy the 
equations 



m 

^^^-i -Q ^ (^■ = l,2,...,m), (19) 



Sjk 



where Sij = l — 2Azj+ZiZj, A = (a^ + 6^ — c^)/2a6 and m is the number of the 
lines. The same equations appear when we try to obtain the lowest energy 
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of the XXZ quantum spin chain. If 1 < A, it is known that zj^s are real and 
Amax = + b^. If A < 1, Zj's he on the unit circle in the complex plane. 
When one consider the limit w ^ oo with fixed m/w, Eq. f|T9|) is reduced to 
an Fredholm integral equation of the second kind, and its solution is obtained 
for m/w = 1/2, with some assumptions. An additional structural transition 
of Ajnax exists at A = —1, where an appropriate parameter yu, A = — cos /i, 
changes from real to pure imaginary. 

These changes correspond to the "phase transitions" of the six- vertex 
model, and of the XXZ quantum spin chain. This fact also means that 
the functional relation for the fractal dimension of this fractal set shows 
structural changeat A = 1 and —1 in the limit w ^ oo. 

If one takes e^j'' = for all i, j and k, then the partition function is equal 
to the number of all the admissible nests of lattice paths for the fixed m. The 
fractal dimension s becomes proportional to the entropy of the nests. This 
case, A = 1/2, corresponds to the high temperature limit of the six- vertex 
model. 

The case A = is particularly simple and Amax is expressed as 

f+°° sinh(f + 7)x 

log Amax = -pei + / -Y—dx, 

J -oo 4a; cosh i^x 

where a/c = sin i(| — 7) and h/c = sin |(| +7). This case is the fractal set 
which corresponds to the XY spin chain. This case also concerns the tiling 
problem (the dimer problem) solved by Kasteleyn[17] and Temper ley and 
Fisher [IB]. The solutions of tiling problems are classified by Cohn, Kenyon 
and Propp[IS] and also classified [T2] with the use of the n-equivalence. 
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lattice models 


fractal geometry 


transfer matrix V 


matrix A*^*^ 


n-equivalence 


transitivity 


free energy 


equation for the fractal dimension s 



Table. 1: Correspondence between lattice models and fractal geometry. 



Figurecaptions: 

Fig. 1(a): Six vertices, corresponding lines and associated energies. 
Fig. 1(b): An allowed line configuration and the transfer matrix V. 
Fig. 2: Dots and directed edges which denote contraction functions. 
Fig. 3: The rectangular lattice with w lows and the order of bonds. 
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